Probing the mechanical properties of graphene using a corrugated elastic substrate 
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The exceptional mechanical properties of 
graphene have made it attractive for nano- 
mechanical devices and functional composite ma- 
terials [Ij. Two key aspects of graphene's me- 
chanical behavior are its elastic and adhesive 
properties. These are generally determined in 
separate experiments, and it is moreover typ- 
ically difficult to extract parameters for adhe- 
sion. In addition, the mechanical interplay be- 
tween graphene and other elastic materials has 
not been well studied. Here, we demonstrate a 
technique for studying both the elastic and ad- 
hesive properties of few-layer graphene (FLG) 
by placing it on deformable, micro-corrugated 
substrates. By measuring deformations of the 
composite graphene-substrate structures, and de- 
veloping a related linear elasticity theory, we 
are able to extract information about graphene's 
bending rigidity, adhesion, critical stress for in- 
terlayer sliding, and sample-dependent tension. 
The results are relevant to graphene-based me- 
chanical and electronic devices, and to the use of 
graphene in composite [Ij, flexible [2j, and strain- 
engineered [3j materials. 

The elastic properties of graphene have previously 
been measured using nano-indentation |1] and pressur- 
ization [5] techniques, and Young's modulus E was found 
to be exceptionally high, ^ 1 TPa. Graphene's van der 
Waals adhesion to surfaces has been examined theoret- 
ically [S], and local adhesion to nanoparticles has been 
studied [7]. Substrate interactions, due to surface adhe- 
sion, highly modify graphene's doping [83 9, and carrier 
mobility. In addition, adhesion to sidewalls in suspended 
graphene introduces strain and modifies mechanical be- 
havior [TOlITT]. Here, we explore both elasticity and adhe- 
sion, which are evident in the interaction between micro- 
scale-corrugated elastic substrates and graphene samples 
of thickness ranging from 1 to 17 atomic layers. By using 
an atomic force microscope (AFM) to determine surface 
adhesion and deformations, we find that the FLG can 
fully adhere to the patterned substrate, and that thicker 
samples flatten the corrugated substrate more than thin- 
ner samples do. By developing a linear elasticity theory 
to model the flattening and adhesion as a function of 
layer thickness, we are able to extract estimates of, or 
bounds on, various fundamental and sample-dependent 
properties of the system. 

Sample substrates were prepared by casting a 3 mm 
thick layer of polydimethylsiloxane (PDMS) — which 
cures into a flexible, rubbery material — onto the exposed 
surface of a writable compact disc. This resulted in ap- 
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FIG. 1. a. AFM 3D topographic image of FLG on corru- 
gated PDMS. AFM lieight and phase data are superimposed 
to create color contrast, b. Illustration of FLG-PDMS inter- 
action, showing how FLG adheres to and flattens the PDMS 
corrugations. The coordinates and amplitudes relevant to the 
theoretical model are labeled. 



proximately sinusoidal corrugations on the PDMS, hav- 
ing a wavelength of 1.5 //m and a depth of 200 nm (see 
Fig. la). Graphene was then deposited onto the PDMS 
via mechanical exfoliation |12| . Candidate samples were 
first located using optical microscopy, then imaged on 
an Asylum Research MFP-3D AFM. Figure la shows a 
topographic image of FLG on the PDMS; it is evident 
from the image that the graphene conforms to the corru- 
gations, as illustrated in Fig. lb. 

In order to fit the experimental data to a theoretical 
model, it was necessary to determine (1) the thickness 
of the FLG, (2) the adhesion between the FLG and the 
PDMS, and (3) the height profile of the PDMS-FLG sys- 
tem. The soft, non-standard substrate created difficulty 
in measuring FLG thickness via established AFM and 
Raman techniques. Thus, the thickness was determined 
by using the AFM (in contact mode) to fold the flake 
onto itself, and then using the AFM again to measure the 
resultant FLG-FLG step height. An example of AFM- 
folded FLG is shown in Fig. 2a. The same flake could be 




FIG. 2. a. Micrograph of FLG folded by AFM tip. FLG- 
FLG edges are used for thickness measurements. Dashed ar- 
row points to fold and shows the direction in which AFM 
tip was dragged, b. AFM phase images of FLG on PDMS. 
Main figure: 2-layer graphene (lower left section of FLG flake) 
showing uniform adhesion over plateaus and valleys, indicated 
by homogeneous contrast over them. Adhesion loss, evident 
as bright and dark lines, occurs only at steepest regions be- 
tween plateaus and valleys. Inset: 13-layer graphene showing 
inhomogeneous phase and adhesion; bubble-like structure in- 
dicates where adhesion changes. Scale bar is 2 fim. 



folded at multiple locations to increase accuracy (typical 
accuracy was 1-2 layers), although, because the method 
was destructive, it had to be undertaken once all other 
measurements were completed. 

The degree of adhesion between the FLG and PDMS 
was obtained by measuring the phase of the oscillation 
of the AFM cantilever. This phase is determined by the 
electrostatic properties of the surface; in other words, 
sections having the same adhesion have common elec- 
trostatic properties and thus a common phase. The 
main image of Fig. 2b shows a two-dimensional phase 
map for 2-layer FLG, in which the phase is uniform 
across the sample (except where adhesion is lost at the 
steepest slopes of the corrugation). These data demon- 
strate the near-conformal adhesion between the FLG and 
the PDMS, and are consistent with previous work on 
graphene placed over more shallow depressions [TT] . The 
AFM height data similarly indicate that the FLG ad- 
heres to the corrugations of the PDMS (e.g., see Fig. 3). 
In contrast, the inset to Fig. 2b shows the phase data for 
13-layer FLG, where bubble-like structure appears across 
the sample, showing that the phase is not uniform and, 
hence, that the FLG does not adhere well to the PDMS. 
In general, we found that samples having fewer than ~11 
layers showed full adhesion, whereas thicker samples did 
not fully adhere. The adhesive properties did not seem 
to depend on the size or aspect ratio of the graphene 
samples, only their thickness. 

The most remarkable aspect of the FLG-PDMS sys- 
tem is that interplay between the rigidity of the graphene 
and the shear modulus of the PDMS causes the FLG to 
become corrugated and the PDMS to be flattened. Fig- 
ures 3a and b show image and height measurements for 
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FIG. 3. Image (top) and height measurements (bottom) for 
a. 8-layer and b. 13-layer graphene. Red lines show trajecto- 
ries of scans over graphene, corresponding to red height curves 
(averaged between the dotted lines). Blue lines show scans of 
surrounding PDMS substrate. Scans over PDMS alone are 
taken far from FLG, to provide a baseline height unaffected 
by FLG. Note the scale difference between vertical and hori- 
zontal axes in height vs. distance curves. 



8- and 13-layer FLG on PDMS, respectively. In Fig. 3a 
it is clear that the FLG maintains the basic shape of the 
PDMS corrugations, but pulls the corrugations up in the 
valleys and pushes them down on the plateaus. In con- 
trast. Fig. 3b shows that 13-layer FLG sits on top of the 
PDMS; while the FLG likely strongly deforms the PDMS, 
the amount of flattening is difficult to determine because 
of the lack of adhesion. Figure 4 shows the fractional 
height difference between the FLG-PDMS composite and 
the bare PDMS, plotted against graphene thickness, for 
18 samples (measured on 9 different PDMS substrates); 
it is clear that the amount of flattening increases with 
layer number. 

We now develop a linear elasticity theory that models 
the effect of stress at the interface between the relatively 
flat FLG (ignoring the nanometer-scale intrinsic ripples) 
and the corrugated PDMS. To start, we consider an un- 
deformed PDMS substrate having a corrugated surface: 
h{x) = Hcoskx (see Fig. lb for coordinates). Height 
proflles that are not simple sine waves can be handled via 
the superposition of suitable terms. Graphene adheres to 
and flattens this surface, reducing the corrugation ampli- 
tude to H — AH (as discussed below). The normal stress 
S cos kx required to deform the graphene in this way fol- 
lows from thin plate theory [T3], and is related to the 
deformed height proflle via 



MV^ [{H - AH) cos kx] ^ S cos kx, 



(1) 



where M is the flexural, or bending, rigidity of the FLG, 
so that S = AIk*{H—AH). An equal and opposite stress 
acts on the PDMS substrate, so we next determine how 
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the PDMS responds to this stress. 
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FIG. 4. Data and fits for flattening factor (fractional defor- 
mation of PDMS) vs. number of layers. Symbols show mea- 
sured flattening for 18 FLG samples. Samples with thick- 
nesses > 11 layers are shown with open triangles, as AFM 
height measurements are likely modified by the lack of ad- 
hesion to the substrate. A value of _F = corresponds to 
no flattening of the PDMS, whereas F — 1 corresponds to 
totally flattened PDMS. Error bars are related to the un- 
certainty in layer number and the spatial non-uniformity of 
flattening. Red curve is least squares flt to model, assuming 
zero tension in samples. Blue curve is predicted flattening for 
samples having zero tension (see text). 

We regard the PDMS as being a semi-infinite, 
isotropic, incompressible medium, and describe distor- 
tions of it by means of a displacement field u{x, z). (We 
only consider configurations that are translationally in- 
variant in the y-direction, and hence are effectively two- 
dimensional.) At the linear level, to which we are work- 
ing, incompressibility implies divergence-free displace- 
ments: V-M = 0. The displacement u obeys the condition 
of mechanical equilibrium, i.e.. 



(2) 



in which p is the pressure field, introduced to implement 
incompressibility, which requires that V^p = 0. In deter- 
mining u, the appropriate boundary condition involves 
specifying the normal component of the stress, which 
amounts to the condition (summing over repeated in- 
dices) 

^ Ui [diUj (f) + djUi (f)) - Uj p (f) = fj (f) , (3) 

in which f{f) is the external force acting on the PDMS 
at its surface and n{r) is the unit normal vector point- 
ing outward from the PDMS surface. The force exerted 
by the graphene on the PDMS can be taken to have a 



simple, sinusoidal form: {fx,fz) = (0, — S'cos fccc). It is 
then straightforward to show that the PDMS responds 
by undergoing the position-dependent displacement 



S 



kz sin kx 
(1 + kz) cos kx 



(4) 



In particular, the surface of the PDMS is displaced ac- 
cording to {ux,Uz)\z=o = — (S'/2/ifc)(0, cos /cx), which cre- 
ates a surface profile {H — AH) cos kx. The amplitude is 
diminished from H by an amount 



AH 



S _ Mk*{H - AH) 
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Rearranging gives the flattening factor 

{Mk^/2n) 



AH 

IT 



1 + (A/fc3/2^)' 



(5) 



(6) 



A model of the FLG as a uniform elastic continuum sug- 
gests that M = Dn? , where n is the number of graphene 
layers and D is the bending rigidity of a single graphene 
layer. We assume the thickness of rt-layer graphene to be 
n times the thickness of single-layer graphene. Thus, we 
arrive at the result that F should increase with increasing 
n. 

We now compare the model to the data, to eluci- 
date the FLG's mechanical behavior. Figure 4 shows 
measured values of F vs. n, along with a least squares 
fit to equation (6) (red curve). From the fit we ex- 
tract a dimensionless graphene rigidity parameter G = 
Dk^/2fi = 0.00091. The shear modulus fi of PDMS was 
measured separately via nanoindentation and an ultra- 
sonic probe, which gave = 0.4 MPa and 0.23 MPa, re- 
spectively (the difference is possibly due to probing sur- 
face vs. bulk moduli). Using /z = 0.4 MPa, as well as 
k — 4.2 /im^^, the bending rigidity of graphene is then 
obtained as £> = 9.8xl0"^*Nm. This value is higher 
than that predicted using Kirchhoff plate theory, from 
which D = Et^/U (l - v"^) = 2.9xlO-i*Nm, using the 
graphene Young's modulus E 900 GPa [13] , Poisson 
ratio 1/ w 0, and single-layer thickness t = 0.335 nm |14j . 
The predicted values for F are plotted as the blue curve 
in Fig. 3. 

The spread in the data is greater than can be accounted 
for by measurement uncertainty. This leads us to hy- 
pothesize that the discrepancy between extracted and 
predicted values of F (or D) is caused by tension in the 
graphene, resulting from sample-dependent friction be- 
tween the PDMS and the graphene. A tension T could 
modify the flattening factor in equation (6), giving 



F . 



AH _ {Mk^/2y) + (rfc/2^) 
l/o~ ^ 1 + (Mfc3/2/i) + (Tfc/2^)' 



(7) 



If we assume that the difference between the predicted 
values of F (blue curve) and the data in Fig. 4 is due to 
tension, we can use equation (7) to extract a value of ten- 
sion for each sample. This yields tension values between 
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and 0.20 Nm""'^, with no discernible trend with thick- 
ness. The tension is positive for each sample, to within 
the margin of error, consistent with the picture that fric- 
tion opposes the contraction of FLG as it conforms to 
the corrugated PDMS. 

The maximum tension corresponds to a maximum ax- 
ial strain of T/nhE = 7.8x 10~^. We can also use the ten- 
sion to estimate the magnitude of the stress due friction: 
assuming the friction acts over a distance d > 10 /im, we 
find a stress T/d of less than 2x10^ Pa. The condition 
that tension be positive, taken together with our data, 
implies that graphene's bending rigidity is no greater 
than (1.6 ± 0.8) x 10""'^* Nm, consistent with predicted 
values. If, as seems reasonable, the tension is negligi- 
bly small for at least one sample, then this bound would 
become an estimate of the graphene rigidity. 

The data shown in Fig. 4 can also be used to estimate 
the normal interface stress S — 2iJ,kFH, which ranges 
from (1.1 to 3.0) x 10^ Pa. The data also show that no 
samples which adhere to the surface have F > 0.6, im- 
plying that the adhesive strength between the graphene 
and the PDMS is < 3.0 x 10^ Pa. Note that this stress is 
model-independent . 

We can extract bounds on the graphene-PDMS ad- 
hesion energy by considering that the energy of the 
adhesion must be at least as large as the spatially- 
averaged elastic deformation energy. This energy can 
be regarded as a sum of contributions: the elastic de- 
formation of the substrate {l/2)aik{uik — {p/^)5ik) — 
SAHcos^{kx)/2; the FLG deformation under tension 
{T/2){Fh'{x)f = T{FHfk'^cos^{kx)/2; and the FLG 

bending {M /2){Fh" {x)f = M{FH)^k'^ cos^{kx)/2 [13]. 
Ignoring the negligible tension contribution, these sum to 
jikFH^ /2. The maximum elastic energy (which is also 
the lower bound of the adhesive energy) in our samples is 
thus 0.044 eV/nm . This is consistent with the theoreti- 
cal prediction of 0.04 eV/nm for the van der Waals adhe- 
sive energy between graphene and a Si02 substrate [T5] . 
Absent any other significant contributions to the energy 
budget (such as work done against friction), so adhesion 
energy must equal the elastic energy and 0.044 eV/nm^ 



becomes an estimate of the adhesive energy. 

Because interlayer forces within FLG are weak, it is 
also important to consider shear forces, which could cause 
the graphene layers to slide relative to one another. In 
this case, the impact on the flattening factor is to replace 
the bending rigidity of a cohesive sample Dir' by the sum 
Tlia ^"-a of bending rigidities Dn^ of the individual 
slabs. We find that such a model does not improve the fits 
to the data, and thus find no evidence that the graphene 
layers slide. The physical effect of sliding would be to 
decrease the flattening factor, and thus sliding cannot ex- 
plain the discrepancy between the theoretical values and 
data of Fig. 4. We hypothesize that there does, however, 
exist some critical shear stress beyond which layers slide 
relative to one another. Considering the finite thickness 
of the FLG, Mindlin plate theory [THj shows that the 
boundary stress needed to deform the FLG generates a 
central shear strain e of 

e^n^t^k^{H - AH)/A. (8) 

The absence of evidence for sliding in samples of < 11 
layers thus suggests a lower bound on the critical shear 
strain ecrit of > 1.2 x 10~^. Multiplying by the graphene 
shear modulus, which we take to be half its Young's mod- 
ulus, this implies a critical shear stress between the layers 
of > 5.6 MPa. 

To conclude, we have developed a method of mea- 
suring the mechanical properties of graphene using de- 
formable, micro-corrugated substrates. We are able to 
put bounds on — or extract estimates for — fundamental 
properties such as graphene's bending rigidity, critical 
shear stress, and the FLG-PDMS adhesive strength and 
energy. We also extracted sample-dependent properties 
such as the tension and normal interface stress. The ex- 
perimental and theoretical techniques developed in this 
paper may be readily extended to various substrates hav- 
ing a range of surface geometries. 
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